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We investigate the electron-transport properties of ethyne-bridged diphenyl zinc-porphyrin molecules sus¬ 
pended between gold (111) electrodes by first-principles calculations within the framework of density func¬ 
tional theory. It is found that the conductance of a molecular junction in which phenyl and porphyrin rings 
are perpendicular is reduced by three orders of magnitude compared with that of a junction in which the 
phenyl and porphyrin rings are coplanar. In the coplanar configuration, electrons are transmitted through 
n states, which extend over the whole molecule. In the perpendicular configuration, the conductance is 
suppressed because of the reduction of electron hopping between n states of the phenyl ring and cr states 
of the porphyrin ring. 


1. Introduction 

Controlling the electron transport through molecules that bridge two electrodes is a basic step 
toward the development of sophisticated nanoscale devices. 1,2) In recent years, the results of 
experimental 3- ^ and theoretical studies 6- ^ have predicted a promising future for nanodevices 
utilizing molecular junctions. Among many of the applications, several schemes have been 
proposed for the design and construction of molecular switches. 5,8,9) The key idea is to find 
two or more distinct states of molecules with significantly different conductance by applying 
an external bias or using a scanning tunneling microscope (STM) tip to manipulate the sys¬ 
tem. 5, 10) In addition, with the advance of laser techniques, intense laser pulses can be used to 
manipulate the torsional motion of some conjugated molecules, which may produce ultrafast 
molecular switches. 11-13) 

Among the classes of conjugated structures in which single molecule charge transport 
measurements have been performed, porphyrin-based molecular structures show exceptional 
electronic structural characteristics. By using the STM break junction method, Li et al. 
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fabricated a molecular junction formed by ethyne-bridged diphenyl zinc-porphyrin (edzp) 
molecules suspended between a gold STM tip and a gold (111) substrate. 14) These systems 
are of interest since it has been reported that the relative angle between the phenyl and por¬ 
phyrin rings could substantially affect the conductance of the molecular junctions. 15) Previous 
first-principles calculations on ethyne-bridged diphenyl porphyrin conjugates demonstrated 
that the coplanar conformation of phenyl and porphyrin rings resulted in a much larger cur¬ 
rent than the perpendicular conformation, and a much current ratio of the on/off states was 
observed. 15) However, the authors focused on free-base porphyrins, whereas in actual exper¬ 
iments, zinc-porphyrins are mainly employed as building blocks to tailar the properties of 
materials. 14,16_20) In addition, since the supercell sizes employed in their calculations were 
small and the molecular geometries were not optimized when the materials were sandwiched 
between gold (111) electrodes, a more accurate study needs to be carried out. To the best of 
our knowledge, a comparison of the transport properties of an edzp molecule with various 
torsional angles between the phenyl and porphyrin rings has not been made. 

In this paper, we examine the electron-transport properties of edzp molecules suspended 
between gold (111) electrodes by first-principles calculations. Our results reveal that the con¬ 
ductance of a molecular junction in which phenyl and porphyrin rings are perpendicular is 
greatly suppressed compared with that of a junction in which the phenyl and porphyrin rings 
are coplanar. The conducting channel in the coplanar conformation is related to n states, 
which extend through the whole molecule. In contrast, the transmission in the perpendicular 
conformation is damped because of negligile coupling between n states of the phenyl ring 
and cr states of the porphyrin ring. The rest of this paper is organized as follows. In Sect. 2, 
details of the computational methods and the models of the molecular junctions are given. 
In Sect. 3, results and discussion of the electron-transport properties for these systems are 
presented. We summarize our findings in Sect. 4. 

2. Computational procedures 

All calculations are performed within the framework of density functional theory using a real- 
space finite-difference approach, 21-25) which makes it possible to determine the self-consistent 
electronic ground state with a high degree of accuracy by a time-saving double-grid tech¬ 
nique. 22-2 ^ The real-space calculations eliminate the serious drawbacks of the conventional 
plane-wave approach, such as its inability to treat electron-transport problems accurately. 

Figure 1 shows the computational models of scattering regions connected to gold (111) 
electrodes. Such two-probe systems are divided into three parts: the left and right electrode 
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regions, which contain three gold (111) layers with 32 atoms per each layer, and the scat¬ 
tering region consisting of the edzp molecule chemisorbed at the fee sites of two gold (111) 
surfaces. 15) We study two representative models, in which the phenyl and porphyrin rings of 
the edzp molecule are coplanar and perpendicular, which will be referred to as the edzpc 
and edzpp models, respectively. Tetragonal supercells are employed with the dimensions 
L x = 19.95 A , L y = 11.52 A, and L z = 35.40 A for the scattering region, where L x and 
L y are the lengths in the x and y directions parallel to the surface, respectively, and L z is 
the length in the z direction. Structural optimizations of the scattering regions are performed 
by relaxing all atoms, except for those in the bottommost layers, until the remaining forces 
are less than 0.05 eV/A, while imposing periodic boundary conditions in all directions. The 
exchange-correlation interaction is treated by the local density approximation 2 ^ based on 
density functional theory, and the projector-augmented wave method 27) is used to describe 
the electron-ion interaction. 

The electron-transport properties of the edzp molecules are then computed using scat¬ 
tering wave functions continuing from one electrode to the other, which are obtained by 
the overbridging boundary-matching method. 22,28) In the calculations of the scattering wave 
functions, the norm-conserving pseudopotentials 29) generated by the scheme proposed by 
Troullier and Martins 30) are employed instead of the projector-augmented wave method. To 
determine the Kohn-Sham effective potential, a supercell is used under a periodic boundary 
condition in all directions, and then the scattering wave functions are computed under the 
semi-infinite boundary condition obtained non-self-consistently. It has been reported that this 
procedure is just as accurate in the linear response regime but significantly more efficient than 
performing computations self-consistently on a scattering-wave basis. 31) The grid spacing is 
set at 0.20 A and the Brillouin zone for the scattering region is sampled at the T point to set 
up the Kohn-Sham effective potential. The conductance at zero temperature and zero bias is 
described by the Landauer-Biittiker formula, 32) G = Go 2,, / jt//j 2 v,7v / , where t (/ - is the trans¬ 
mission coefficient at which electrons are transmitted from an initial mode j to a final mode i 
inside the scattering region, v,- and vj are the longitudinal components of the group velocity in 
modes i and j, and Go = e 2 /h, with e and h being the electron charge and Planck’s constant, 
respectively. 

3. Results and discussion 

The calculated conductances of the edzpc and edzpp models, as functions of the energy of 
incident electrons from the left eletrode, are presented in Figs. 2(a) and 2(b), respectively. A 
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feature of immediate interest is that the conductance is suppressed by three orders of mag¬ 
nitude when the phenyl and porphyrin rings are perpendicular to each other. At the energy 
of E F - 0.35 eV, where E F is the Fermi level, the conductance of the edzpc model has one 
broad peak with a height of 1.00 Go while that of the edzpp model has one narrow peak with 
a height of 3.48 x 10 _3 G 0 . 

In Fig. 3, we plot distribution of the local density of states (LDOS) of the two models, 
which are plotted by integrating them along the x-y plane, p(z,E ) = \ifr(r, E)\ 2 dr\\, where 
r = (x, y, z), is the wave function, and E is the energy of the states. In the energy window 
from E F - 0.50 eV to Ep, there are electronic states localized at the porphyrin ring, which 
correspond to the observed peaks in the conductances of the two models. From E F to Ep+0.50 
eV, the conductances are almost zero owing to the presence of the gap in the LDOS. 

It is notable that the conductance of the edzpp model is very small although there is a large 
density of states at the porphyrin ring from Ep- 0.50 eV to Ep. To understand the effect of the 
rotation of the phenyl rings on the electron-transport property, in Fig. 4 we plot the projected 
density of states (PDOS) on the p x and p y atomic orbitals of the carbon atoms indicated 
by the arrows in Fig. 1. Around the Fermi level, n states are expected to contribute to the 
electron transport. 15) In the edzpc model, the electrons are propagated from the left electrode 
by states with the p x character. Since there is a substantial overlap between the PDOS on the 
p x orbital of the A cop atom of the phenyl ring and that of the B cop atom of the porphyrin ring, 
the conductance is large. In contrast, when the phenyl rings are rotated by 90°, the electrons 
are propagated from the left electrode by states with the p y character. Because the overlap 
between the PDOS on the p y orbital of the A per and B per atoms is small, the electron transport 
is reduced. 

Finally, to investigate which states contribute to the electron transport in more detail in 
Fig. 5, we depict the isosurfaces of the charge density distribution of the scattering waves for 
electrons injected from the left electrode. The peak in the conductance of the edzpc model 
corresponds to the transmission of n states of the molecule. Because of the orthorgonality 
between the phenyl and porphyrin rings in the edzpp model, n states of the porphyrin ring do 
not conduct electrons whereas cr states of the porphyrin ring can contribute to the electron 
transport. However, electrons are reflected at the junctions between the phenyl and porphyrin 
rings owing to the small electron-hopping probability from n states of the phenyl ring to cr 
states of the porphyrin ring, which results in the small conductance. 
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4. Conclusion 

We have studied the electron-transport properties of edzp molecules suspended between gold 
(111) electrodes by first-principles calculations. It is found that the conductance through the 
molecular junction is greatly suppressed when the phenyl rings are perpendicular to the por¬ 
phyrin ring. The peaks of conductances of 1.00 Go and 3.48 x 10 -3 Go for the edzpc and 
edzpp models, respectively, are observed at the energy of E F - 0.35 eV. The charge density 
distributions of the scattering waves reveal that the conduction through the edzpc model origi¬ 
nates from the transmission of n states of the porphyrin ring while the conduction through the 
edzpp model is related to <j states of the porphyrin ring because n states of the porphyrin ring 
do not contribute to the electron transport. Our results suggest that this molecular junction 
can be used to construct molecular switches. 
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Fig. 1 . (Color) Computational models of (a) edzpc and (b) edzpp molecules. Red, green, blue, white, 
purple, and gold spheres indicate zinc, carbon, nitrogen, hydrogen, sulfur, and gold atoms, respectively. 
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Fig. 2. Conductances of (a) edzpc and (b) edzpp models as functions of energy of incident electrons. Zero 
energy is chosen to be at the Fermi level. 
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Fig. 3. (Color) Distributions of LDOS of (a) edzpc model and (b) edzpp model integrated on plane parallel 
to electrode surface as functions of relative energy from the Fermi level. Zero energy is chosen to be at the 
Fermi level. Each contour represents twice or half the density of the adjacent contour lines, and the lowest 
contour is 2.78 x 10 -5 e/eV/A. The atomic configurations are given as a visual guide below the graph, and the 
symbols have the same meanings as those in Fig. 1. 
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Fig. 4. (Color online) PDOS of p x and p y orbitals of (a) A cop and (b) B cop atoms of edzpc model shown in 
Fig. 1(a), and PDOS of (c) A per and (d) B per atoms of edzpp model shown in Fig. 1(b). Zero energy is chosen 
to be at the Fermi level. 
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Fig. 5. (Color) Top and side views of isosurfaces of charge density distributions of scattering waves emitted 
from left electrode of (a) edzpc model and (b) edzpp model at the energy of E F - 0.35 eV. The iso value is 
1.24 x 10 -3 e/A 3 /eV. The symbols have the same meanings as those in Fig. 1. 












